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Abstract 

We study the cosmological evolutions of the equation of state for dark energy wde in the expo- 
nential and logarithmic as 'well as their combination /(T) theories. We sho'w that the crossing of 
the phantom divide line of wde = — 1 can be realized in the combined /(T) theory even though 
it cannot be in the exponential or logarithmic f{T) theory. In particular, the crossing is from 
''^DE > — 1 to uiDE < —1) in the opposite manner from f{R) gravity models. We also demonstrate 
that this feature is favored by the recent observational data. 
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(CMB) radiation 
and weak lensing 



I. INTRODUCTION 

Cosmic observations from Supernovae la (SNe la) cosmic microwave background 
2-4|, large scale structure (LSS) 5|, baryon acoustic oscillations (BAO) jo], 
7| have implied that the expansion of the universe is currently accelerating. 
This is one of the most important issues in modern physics. Approaches to account for the 
late time cosmic acceleration fall into two representative categories: One is to introduce 
"dark energy" in the right-hand side of the Einstein equation in the framework of general 
relativity (for a review on dark energy, see 8|). The other is to modify the left-hand side of 
the Einstein equation, called as a modified gravitational theory, e.g., f{R) gravity joHl. 

As another possible way to examine gravity beyond general relativity, one could use 
the Weitzenbock connection, which has no curvature but torsion, rather than the curvature 
defined by the Levi-Civita connection. Such an approach is referred to "teleparallelism" (see. 



e.g., |12|-|l5|). which was also taken by Einstein 16|. To explain the late time acceleration of 
the universe, the teleparallel Lagrangian density described by the torsion scalar T has been 



extended to a function of T 
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ISp. This idea is equivalent to the concept of f{R) gravity. 



in which the Ricci scalar R in the Einstein-Hilbert action is promoted to a function of R. 
Recently, f{T) gravity has been extensively studied in the literature 
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25|. 



In this paper, we explicitly examine the cosmological evolution in the exponential /(T) 



theory 18|, |23| in more detail with the analysis method in Ref. 26|. In particular, we study 



the equation of state {vjjm) and energy density (pde) for dark energy. The recent cosmologi- 



cal observational data [27| seems to imply a dynamical dark energy of equation of state with 
the crossing of the phantom divide line Wde = — 1 from the non-phantom phase to phantom 
phase as the redshift z decreases in the near past. However, we illustrate that the universe 
with the exponential /(T) theory always stays in the non-phantom (quintessence) phase or 
the phantom one, and hence the crossing of the phantom divide cannot be realized j23|. It 
is interesting to mention that such an exponential type as f{R) gravity models has been 
investigated in Refs. 28l-l30| . We also present a logarithmic f{T) theory and show that it 
has a similar feature as the exponential one. Our motivation in this paper is to build up 
a realistic f(T) theory in which the same behavior of the crossing of the phantom divide 



^ Models based on modified teleparallel gravity for inflation have been investigated in Ref. [3] • 
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as indicated by the data can be achieved. For this purpose, we will construct a combined 
/(T) theory with both logarithmic and exponential terms. Furthermore, we examine the 
observational constraints on the combined f{T) theory by using the recent observational 
data of SNe la, BAO and CMB. We note that two /(T) models with the crossing of the 



phantom divide have also been proposed in Ref. [22|. 

The paper is organized as follows. In Sec. II, we introduce /(T) theory and derive the 
gravitational field equations. In Sec. Ill, we investigate the cosmological evolutions in the 
exponential and logarithmic f{T) theories. In Sec. IV, we propose a model which combines 
the two theories in Sec. Ill to achieve the crossing of the phantom divide. In Sec. V, we 
explore the observational constraints on the combined model in Sec. IV. Finally, conclusions 
are given in Sec. VI. 

II. THE COSMOLOGICAL FORMULAE IN f(T) GRAVITY 

In the teleparallelism, orthonormal tetrad components e^lx'^) are used, where an index 
A runs over 0, 1, 2, 3 for the tangent space at each point of the manifold. Their relation 
to the metric g'^'^ is given by 

g^. = VABeie^, (2.1) 

where fi and u are coordinate indices on the manifold and also run over 0,1,2,3, and 
forms the tangent vector of the manifold. 

The torsion T^^^, and contorsion K'^'^ p tensors are defined by 

TV = {d,et - d^e^) , (2.2) 
K''^ = {T^^p - T^^p - Tp . (2.3) 

Instead of the Ricci scalar R for the Lagrangian density in general relativity, the teleparallel 
Lagrangian density is described by the torsion scalar T, defined as 



where 



(2.4) 



S^^l {K'\ + - T^'a) ■ (2.5) 



Consequently, the modified teleparallel action for /(T) theory is given by 18 1 



d'x\e\iT + fiT)) , (2.6) 



where |e| = det (e^) = ^J—g. 

We assume the flat Friedmann-Lemaitre- Robertson- Walker (FLRW) space-time with the 
metric, 

dg" = dt^ - a^{t)dx^ , (2.7) 

where a{t) is the scale factor. In this space-time, g^y = diag(l, — a^, — a^, — a^) and therefore 
the tetrad components = (l,a,a, a) yield the exact value of torsion scalar T = —6H^, 
where H = d/a is the Hubble parameter. We use units of = c = h = 1 and the 
gravitational constant G = Mp^ with the Planck mass of Mpi = 1.2 x 10^^ GeV. 

It follows from the variation principle that in the flat FLRW background, the modified 
Friedmann equations are given by [3] 

= 2H'fT , (2.8) 

, _ WttGPm + 6H' + f + UH'fr 

^ > ~ 2AH^fTT-2-2fT ' ^ ^ 

where a prime denotes a derivative with respect to In a, fr = df(T)/dT, Jtt = d'^fiT)/dT^, 
and pm and Pm are the energy density and pressure of all perfect fluids of generic matter, 
respectively. 

By comparing the above modified Friedmann equations (12.81) and (12. 9p with the ordinary 
ones in general relativity: 

H' = ^ipM + PDE) , (2.10) 
{H^y = -StiG (pm + Pm + Pde + Pde) , (2.11) 

the energy density and pressure of the effective dark energy can be described by 

PDE = (-/ + 2T/^) , (2.12) 

P - ^ / - + 2T^fTT . . 

^^'^-levrG l + fT + 2TfTT ■ ^ ' 

The equation of state for dark energy is defined as 1^ 

_PdE_ T' fT + 2TfTT _ f/T-fT + 2TfTT , 

"^""^ PDE ^3Tf/T-2fT (1 + /t + 2T/tt)(//T-2/t)- ' 

Since we are interested in the late time universe, we consider only non-relativistic matter 
(cold dark matter and baryon) with pu = Pm and Pm = -Pm = 0, where pm and Pm are 
the energy density and pressure of non-relativistic matter, respectively. Consequently, from 



Eqs. (12. 8 p and (12. 9p it can be shown that the effective dark energy satisfies the continuity 
equation 

= Pde = -3(1 + wde)pde, (2.15) 

where = In a. 



III. COSMOLOGICAL EVOLUTION IN THE EXPONENTIAL AND LOGA- 
RITHMIC /(T) THEORIES 

To analyze the cosmological evolution of the equation of state for dark energy wj^e in 



/(T) theory, we define a dimensionless variable 26 1 



yH = — -a 3.1 



where 



m = , (3.2) 

and pi? = Pm{z = 0) is the current density parameter of non-relativistic matter with the 
redshift z = 1/a — 1. Using Eq. (I2.15p . we obtain the evolution equation of the universe as 
follows: 

y'^ = -S(l + WBE)yH ■ (3.3) 

Note that wbe is a function of T, while T is a function of H^. Moreover, it follows from 
Eq. dSlD that = rh^ [yu + a"^). 

A. Exponential /(T) theory 



We consider the exponential f{T) theory in Ref. [18|, given by 

/(T) = aT (1 - eP^°/^) (3.4) 

with 

1 

a = ; r — , (3.5) 

1 - (1 - 2p) eP ^ ^ 

where p is a constant with p = corresponding to the ACDM model and Tq = T{z = 0) 
is the current torsion. Here, f2m^ = PmVPcrit "where is the energy density of non- 
relativistic matter at the present time and p^^l^ = SH^/ (SttG) is the critical density with 
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Wde 



Wde 




FIG. 1: Wde as a function of the redshift z for \p\ =0.1 (solid line), 0.01 (dashed line), 0.001 
(dash-dotted line) and = 0.26 in the exponential /(T) theory, where the left and right panels 
are for p > and p < 0, respectively. 

Hq being the current Hubble parameter. We note that a in Eq. (13. 5p has been derived from 
^ pg^/p(2t = = -a [1 - (1 - 2p) eP] by using Eq. (12321) with T = Tq and 

Pde ~ Pde{z = 0). We remark that the theory in Eq. (13.41) contains only one parameter 
p if the value of r2m^ is given. The values of other dimensionless quantities at z = can 
be calculated by using p and ^Im , e.g., m^/To = ^SvrGpmVsj / {—^Hq) = —0,^^/6 and 
yniz = 0) = Hl/m^ - 1 = l/fi^^^ - 1. 



In Fig. [H we depict the equation of state for dark energy wde in Eq- ( 2.14p as a function 
of the redshift z for \p\ = 0.1, 0.01, 0.001. We have used fim'' = 0.26 4| in all numerical 
calculations. From Fig. [H we see that wpy does not cross the phantom divide line wbe = 
— 1 in the exponential /(T) theory 18|, |23|]. In particular, for p < the universe always 
stays in the non-phantom (quintessence) phase (wde > —1), whereas for p > it in the 
phantom phase {wbe < —1)- The present values of wbe are wbe{z = 0) = —1.03, —1.003, 
-1.0003, -0.954, -0.996 and -0.999 for p = 0.1, 0.01, 0.001, p = -0.1, -0.01 and -0.001, 
respectively. The larger \p\ is, the larger the deviation of the exponential /(T) theory from 
the ACDM model is. We note that in solving Eq. (13. 3p numerically, we have taken the initial 
conditions at ^ = as yniz = 0) = 2.8. 
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FIG. 2: Wde as a function of |T/To|, where the thm sohd hne shows wde = —1 (cosmological 
constant). Legend is the same as Fig. [TJ 




FIG. 3: IT/TqI as a function of the redshift z for \p\ = 0.1. 



We display wde as a function of |T/To| in Fig. [2], where the left and right panels are for 
p > and p < 0, respectively. In Fig. [3], we show the cosmological evolution of \T/To\ as a 
function of the redshift z for \p\ =0.1. The qualitative behaviors of \T/Tq\ for \p\ = 0.01 
and 0.001 are similar to those for \p\ = 0.1. In addition, the cosmological evolution of 
p^p. = pde/ Pde ^ function of the redshift z is shown in Fig. HI 

From Fig. [2], we find that the crossing of the phantom divide occurs around \T/To\ ^ 
0.7 < 1. This means that the crossing cannot be realized in the past and far future, unlike 
those in f{R) models [29|. Explicitly, by denoting |T/To|cross being the crossing point, we 
have |T/To|cross ^ 0.772, 0.744, 0.740, 0.689, 0.736, and 0.740 for p = 0.1, 0.01, 0.001, -0.1, 
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FIG. 4: = Pde/Pde ^ function of the redshift z, where the thin sohd hne shows wbe = — 1 
(cosmological constant). Legend is the same as Fig. [TJ 



—0.01, and —0.001, respectively. Clearly, the universe reaches the line wde = —1 without a 
crossing even in the very far future {z ^ —1). 

In the high z regime, the universe is at the matter-dominated stage and therefore, |T| = 
QH^ = IGttG (pm + Pde) decreases with time. Around the present time {z = 0), pm is smaller 
than Pde and dark energy becomes dominant over matter. Fig. [3] illustrates that |T/To| 
decreases monotonously with z. At the infinite future, |T/To|cross is reached. However, the 
the universe is in either the phantom phase [p > 0) or non-phantom phase {p < 0) without 
the crossing. Since by that time the contribution from non-relativistic matter and radiation 
is supposed to be negligible, such an evolution pattern of \T/To\ is driven by the /(T) term. 
In the future {z < 0), Pde arrives at the crossing point and finally becomes constant as 
shown in Fig. |H 

The right-hand side of the first equality in Eq. ( I2.14p implies that when wde crosses the 
phantom divide line of wde = —1, the term T' {fx + 2T/tt) / [3T (//T — 2/t)] has to flip 
the sign. Since the sign of T is fixed to be negative due to T = —6H^, the sufficient and 
necessary conditions for the crossing of the phantom divide is that the sign of the combination 
T' [fr + '^T fxr) / {f/T — 2 fx) changes in the cosmological evolution. Note that the sign of 
H is opposite to T' since T' = —12HH' = —12H. When the energy density of dark energy 
becomes perfectly dominant over those of non-relativistic matter and radiation, one has 
wj)E ~ Wes = — 1 — 2H / (3if^) = Ptot/ptot, where Wes is the effective equation of state for 
the universe {q], and ptot = PDE + Pm + Pr and Ptot = PoE + Pr are the total energy density and 
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pressure of the universe, with and Pj. being the energy density and pressure of radiation, 
respectively. For if < (> 0), Wc?i > —1 (< —1), representing the non-phantom (phantom) 
phase, while Wcs = —1 for H = 0, corresponding to the cosmological constant. Clearly, 
the physical reason why the crossing of the phantom divide appears is that the sign of H 
changes from negative (the non-phantom phase) to positive (the phantom phase) due to the 
dominance of dark energy over non-relativistic matter and radiation. 

We now investigate why the crossing of the phantom divide cannot occur in the expo- 
nential f(T) theory in Eq. (13.41) for p < 1. From Eq. (13. 4p . we obtain 

fT = a(^l- e^^«/^ + ^e^^«/^^ , (3.6) 

/TT = -a(f^)'^e^^«/^ (3.7) 

To illustrate our results, we only concentrate on the limits of < p ^ 1 and X = pTo/T <^ 1 
In this case, Tq/T < 1, which corresponds to the region from the far past to the near future. 
Consequently, Eqs. (13. 4p . (13. 6 p and (13. 7p are approximately expressed as 

« X + — , /t^V' TfTT^-aX\ (3.8) 



T V 2 

By using these equations, we find 



l-2fT^-ax(l + ^] , (3.9) 



T V 2 

/t + 2T/tt^-^. (3.10) 

From Eqs. (13.91) and (I3.10p . we see that the signs of both f /T — 2 fx and fr + 2T fxr do not 
change. Note that the sign of T' is also unchanged because \T /Tq\ decreases monotonously 
with z until |r/To|min as shown in Fig. [3l As a result, the crossing of the phantom divide 
cannot be realized in the exponential /(T) theory. 

In Fig. ISl we illustrate the fractional densities of dark energy (^de = P'D'e/ Pcltli non- 
relativistic matter (f2m = Pm/Pcrit) ^n*^ radiation (fij. = Pr/Pcrit) functions of the redshift 
z for \p\ = 0.1. The cosmological evolutions of ^de? and for \p\ = 0.01 and 0.001 are 
similar to those for |p| = 0.1. In order to analyze not only non-relativistic matter and dark 
energy but also radiation, we have included the contribution from radiation in Eq. (13.11) and 



used the following variable 2^^. 



Jj2 

yu = — - - Xa"^ , (3.11) 
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FIG. 5: r^DE (dashed line), (solid line) and Or (dash-dotted line) as functions of the redshift z 
for \p\ =0.1. Legend is the same as Fig. [TJ 



with 



X 



3.1 X 10" 



(3.12) 



where pi^^ is the energy density of radiation at the present time. 

In the high z regime {z > 3225), the universe is at the radiation-dominated stage (fij. 3> 
fl-DE, > ^m)- As z decreases, the universe enters the matter-dominated stage (^2^ > ^de, 
fijn ^ Qj.). After that, eventually dark energy becomes dominant over matter for z < z^-e, 
where zbe is the crossover point in which Qy)E = ^m- Explicitly, we have zbe = 0.40 
and 0.44 for p = 0.1 and —0.1, respectively. At the present time [z = 0), as an initial 
condition we have used = (0.74,0.26,8.1 x 10"^). Note that z^e = 0.42 in 

the ACDM model. As a result, the current accelerated expansion of the universe following 
the radiation-dominated and matter-dominated stages can be achieved in the exponential 
/(T) theory. 

B. Logarithmic /(T) theory 



In this subsection, we examine a logarithmic f(T) theory, given by 

/(^) = /?^„(f)"'^„(f 



with 



/5 



1 - 

2g-i/2 
10 



(3.13) 



(3.14) 



Wde 



Wde 
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FIG. 6: WDE as functions of z (left panel) and IT/TqI (right panel) for q = 1 and f^m^ = 0.26 in 
the logarithmic f{T) theory, where the thin solid line shows tfoE = — 1 (cosmological constant). 



where g is a positive constant. We note that the theory in Eq. (13.131) contains only one 
parameter q if the value of f2m^ is obtained in the same way as the exponential /(T) theory. 

In Fig. [HI we draw w^e as functions of z (left panel) and \T/Tq\ (right panel) for q = 1 
and f2m^ = 0.26 in the logarithmic f{T) theory in Eq. (13.131) . From Fig. El we observe that 
Wde does not cross the phantom divide line wde = — 1, similar to the exponential theory. In 
this logarithmic theory, the universe is always in the non-phantom phase (u'de > ~1) a-^d 
\T/To\ decreases monotonously with z. The present value of wbe is wj)e{z = 0) = —0.79 for 
q = 1. The right figure in Fig. |6] indicates that the crossing of the phantom divide occurs 
at |T/To|cross = 0.547600 and the minimum of \T/Tq\ is |r/To|min = 0.547600. Thus, the 
universe does not cross the line of Wbe = —1 even in the very far future {z ^ —1). 

We remark that wde is qualitatively independent of the value of q. This is because in 
the logarithmic f{T) theory, wde and p^E can be expressed as 

1 



Wbe 



(0) 



and 



T 



Pde 



167rG 



(1 - m 



T 
To 



-1/2 



-1/2 



(3.15) 



(3.16) 



respectively. Similar to the exponential /(T) theory, it is easy to show that the crossing of 
the phantom divide cannot occur in the logarithmic /(T) theory. 
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IV. REALIZING THE CROSSING OF THE PHANTOM DIVIDE IN /(T) THE- 
ORY 



As shown in the previous section, the universe always stays in the phantom and non- 
phantom phases in the exponential {p > 0) and logarithmic /(T) theories, respectively. 
In this section, we investigate the cosmological evolution in a combined /(T) theory with 
both logarithmic and exponential terms and check if the crossing of the phantom divide can 
happen in this combined theory. The explicit form of the theory is given by 



f{T) = 7 

with 



To(f^)"''ln(!^)-T(l-e-o/-) 



(4.1) 



1 - 

^ ~ 2u-V2 + [i_(i_2n)e"] ' ^ ' ' 

where m is a constant. Here, we have taken p = q = u>Oto simplify our discussions when 
we combine Eqs. (13.41) and (13.131) . We note that the model in Eq. (14. ip contains only one 
parameter u if one has the value of Qm^ like the exponential and logarithmic models. 

In Fig. [71 we plot the equation of state for dark energy as a function of the redshift 
z for n = 1, 0.8 and 0.5 and f2m^ = 0.26 in the combined /(T) theory in Eq. (14. ip . As seen 
from Fig. [71 it is clear that Wde can cross Wde = — 1 at 2 = Across- The universe evolves from 
the non-phantom phase {wbe > — 1) to the phantom one {wde < —1) and finally approaches 
Wde = —1- The present values of wde are wde(2 = 0) = —1.1, —1.1 and —0.94 for u = 1, 
0.8 and 0.5, respectively. The crossing values of 2;cross are 0.70 and 0.36 for u = 1 and 0.8, 
respectively. It is interesting to note that such a crossing behavior is consistent with the 
recent cosmological observational data but it is opposite to those in the viable f{R) 
theories 2^. We note that for u < 0.5, the universe asymptotically approaches wbe = — 1 
without crossing it even in the very far future. The cosmological evolution of wde as a 
function of \T/Tq\ is shown in Fig. [HI while that of \T/Tq\ as a function of z for n = 1 is 
depicted in Fig. [9] In addition, p^^ = Pde/Pde ^ function of z is given in Fig. [TUl 



12 



Wde 
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z 



FIG. 7: Wde as a function of the redshift z ior u = I (solid line), 0.8 (dashed line), 0.5 (dash-dotted 
line) and = 0.26 in the combined /(T) theory with both logarithmic and exponential terms 
in Eq. (|4.ip . where the thin solid line shows wde = — 1 (cosmological constant). 

IVde 




0.5 1 1.5 2 2.5 3 3.5 4 

\T/To\ 

FIG. 8: Wde as a function of \T/Tq\. Legend is the same as Fig. [71 



We now demonstrate that the crossing of the phantom divide can occur in the combined 
/(T) theory in Eq. (14. ip . In this theory, both signs of T and T' are unchanged hke the 
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FIG. 9: |r/To| as a function of the redshift z for u = 1. The range of —0.95 < z < —0.40 is 
magnified in the small window. Legend is the same as Fig. [71 




FIG. 10: = Pde/Pde ^ function of the redshift z. Legend is the same as Fig. [71 
exponential /(T) theory. From Eq. (14. ip . we obtain 

= -2^ (i - 1 ,„ (!^^) ) - , (i - e--n + . 

/ ■T' \ 3/2 / rri \ / rri \ 3 
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FIG. 11: r^DE (dashed line), Q.„i (solid line) and (dash-dotted line) as functions of the redshift 
z for u = 1 . Legend is the same as Fig. [71 



leading to 



2 T 



/t + 2T/TT ~ -7 



T 



T 

1 3 /nTo 
n ~ 2 V ^ 



3/2' 



(4.4) 



where we have only concentrated on the case of uTq/T < 1 to simplify our discussion. From 
Eq. (14. 4p . we see that the sign of f /T — 2 fx is fixed, while that of fr + 2TfTT can change 
when uTq/T becomes larger or smaller than the value of uTq/T = [2/ (3n)]^^^. Hence, the 
sufficient and necessary conditions for the crossing of the phantom divide in the combined 
/(T) theory can be satisfied. 



In Fig. [m we demonstrate the fractional densities of dark energy (^de = PoE/Pcnt)' 
non-relativistic matter (f2m = Pm/Pcrit) radiation [Qj- = Pi/p^^lf) as functions of z for 
u = 1. The evolutions of \T/Tq\, Qbe, and Qj- for u = 0.8 and u = 0.5 are similar to those 
for u = 1. In the high z regime {z > 3225), the universe is at the radiation-dominated stage 
{Qj. ^ ^DE, > ^m)- As z decreases, the matter-dominated stage {fl^ > ^be, ^ ^r) 
follows, and eventually dark energy becomes dominant over matter for z < zde- Explicitly, 
we have z^e = 0.36, 0.40 and 0.49 for u = 1, 0.8 and 0.5, respectively. As a consequence, 



15 



after the radiation and matter-dominated stages, the current accelerated expansion of the 
universe can be reahzed in the combined /(T) theory in Eq. (14. ip . similar to that in the 
exponential f{T) theory. 

Finally, we remark that since the gravitational field equation in f{R) gravity has the 
fourth-order nature, it is possible to reproduce an arbitrary background expansion history. 



provided that we do not consider the stability issue [3l|, |32|). On the other hand, in the 
modified Gauss-Bonnet or f{G) gravities, the ACDM background expansion history cannot 
be achieved all the way down to the far future 33|]. In f{T) gravity, there exists a similar 
restriction. Physically, in such modified gravity theories the modified background expansion 
rate H is in general a function of H itself and H, which means that arbitrary expansion 
histories are not guaranteed to be realizable. 



V. OBSERVATIONAL CONSTRAINTS ON THE COMBINED /(T) THEORY 

In Sec. IV, we have constructed the combined f(T) theory with both logarithmic and 
exponential terms to realize the crossing of the phantom divide. Since it is difficult to 
obtain the information on the form of /(T) from more fundamental theories, it is important 
to examine whether this theory is compatible with the recent observational data. In this 
section we demonstrate how well the combined /(T) theory can fit the observational data. 
We examine the constraints on the model parameter u and the current fractional non- 
relativistic matter density Qm^ by taking the method for the recent observational data. 
We use type la supernovae (SNe la) data from the Supernova Cosmology Project (SCP) 
Union2 compilation [3J|, the baryon acoustic oscillations (BAO) data from the Two-Degree 
Field Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey Data Release 
7 (SDSS DR7) 35], and the cosmic microwave background (CMB) radiation data from 
Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) observations J]. The detailed 



analysis method 
Appendix A. 



30| for the observational data of SNe la, BAO and CMB is summarized in 



We plot the contours of the 68.27% (la), 95.45% (2cr) and 99.73% (3a) confidence lev- 
els (CL) in the [^lm\u) plane from the SNe la, BAO and CMB data for the combined 
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FIG. 12: Contours of 68.27% (la), 95.45% (2a) and 99.73% (3a) confidence levels in the {Q^^\u) 
plane from SNe la, BAO and CMB data for the combined f{T) theory, where the plus sign depicts 
the best-fit point. 

TABLE I: The best-fit values of u, h and Xmin ^^r the combined /(T) and ACDM models. 



Model 


U n^^^ h Xmin 


f{T) 


0.829 0.282 0.691 544.56 


ACDM 


0.275 0.707 545.23 



lis figure, we have fitted three parameters of u, f2m^ and 



/(T) theory in Fig. [H In t 

h = ifo/100/[kmsec~^ Mpc~^] [36|]. The model parameter u is constrained to 0.6 < u < 1.13 
(95.45% CL) whereas 0.255 < Q^^'' < 0.312 (95.45% CL), which is consistent with the 
current observations. The best-fit value, i.e., the minimum x^, in the parameter space is 
xLn = 544.56 with u = 0.829, fi^ = 0.282 and h = 0.691. In Table [D we present the 
best-fit values of u, Qm^ and h, and Xmm combined /(T) theory along with those for 

the ACDM model. From Tabled we see that Xmm of combined /(T) theory is slightly 
smaller than that of the ACDM model. This implies that the combined f(T) theory can fit 
the observational data well. 

We remark that for around the best-fit value of the model parameter u ~ 0.8, wde can 
cross the phantom divide line of wbe = — 1 from the non-phantom phase (wde > —1) to the 
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phantom one (wde < —1) in the combined f{T) theory as shown in Fig. [TJ and therefore the 
observational data may favor the crossing of the phantom divide. This consequence agrees 
with the one in Ref. 22[ |. 



VI. CONCLUSIONS 

We have investigated the cosmological evolution in the exponential f{T) theory, which is 
one of the simplest f(T) theories as it contains only one model parameter p along with the 
current fractional density of non-relativistic matter fim^ • We have explicitly shown that the 
phase of the universe depends on the sign of the parameter p, i.e., for p < 0(> 0) the universe 
is always in the non-phantom (phantom) phase without the crossing of the phantom divide. 
We have presented another simplest /(T) model, which is the logarithmic type. Similar to 
the exponential model, the logarithmic one has only one free parameter q, and it does not 
allow the crossing of the phantom divide. To realize the crossing of the phantom divide, 
we have constructed a f(T) theory by combining the logarithmic and exponential terms. In 
particular, we have shown that the crossing in the combined f{T) theory is from Wbe > — 1 
to Wbe < —1, which is opposite to the typical manner in f{R) gravity models. Furthermore, 
we have also illustrated that this combined theory is consistent with the recent observational 
data of SNe la, BAO and CMB. 
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Appendix A: Analysis method for the observational data 



In this appendix, we explain the analysis method 30| for the observational data of type la 



supernovae (SNe la), baryon acoustic oscillations (BAO) and cosmic microwave background 



(CMB) radiation (for more detailed explanations on the data analysis, see, e.g., |37]). 



1. Type la Supernovae (SNe la) 

The information on the luminosity distance as a function of the redshift z is given by 
SNe la observations. The theoretical distance modulus /ith is defined by 

IJ'thizi) = 5 logio DL{z^i) + /io , (Al) 



where /iq = 42.38 - Slog^o/i with h = /Jo/100/[kmsec"^ Mpc"^] |36|. The Hubble-free 
luminosity distance for the flat universe is described as 



dz' 

Dl{z) = {1 + z) I —-, (A2) 







where E{z) = H{z)/Hq with 



H{z) = Ho\Jnt^ (1 + zf + nf'^ (1 + zf + (1 + ^)3(i+"'de) _ (A3) 

Here, fir^^ = ^^7°^ (1 + 0.2271A''cfr), where fi^"'' is the present fractional photon energy density 
and A^'cff = 3.04 is the effective number of neutrino species {4]. We note that H[z) is evaluated 
by using numerical solutions of Eq. ( 13. 3p . 
The of the SNe la data is given by 

XsN — 2^ 2 ' y^^) 

where /lobs is the observed value of the distance modulus. In what follows, subscriptions 
"th" and "obs" denote the theoretically predicted and observed values, respectively. Xsn 
should be minimized with respect to /iq, which relates to the absolute magnitude, because 



the absolute magnitude of SNe la is not known. i^i Eq. flA4l) is expanded to be [38 1 



xl^ = A- 2fioB + fxlC , (A5) 
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with 



The minimum of Xsn with respect to /xq is expressed as 



?1 

c 



(A6) 



(AT) 



In our analysis, we take Eq. ( ]A7p for the minimization and use the Supernova Cosmology 
Project (SCP) Union2 compilation, which contains 557 supernovae 3J], ranging from z = 
0.015 to z = 1.4. 



2. Baryon Acoustic Oscillations (BAO) 



The distance ratio of = rs{zd)/ Dv{z) is measured by the observation of BAO. Here, 
Dy is the volume- averaged distance, is the comoving sound horizon and z^ is the redshift 
at the drag epoch 35|. The volume-averaged distance Dv{z) is defined as j6| 

nl/3 



Dv(z) 



:i + z?Dl(z ' 



Hz) 



(A8) 



where Da{z) is the proper angular diameter distance for the fiat universe, defined by 

1 r dz' 



Da(z) 



(A9) 



The comoving sound horizon rs{z) is given by 

-1/(1+^) 



da 



am{z' = 1/a - + (3fir/4^^f ) « 



(AlO) 



where fi^"-* = 2.2765 x 10 ^/i ^ and ^l^-y'' = 2.469 x 10 ^/i ^ are the current values of baryon 



)(0) 



51,-2 



and photon density parameters, respectively 

1291(fiL°^/l2)0-251 



J] . The fitting formula for z^j^ is given by 



1 + 0.659(^]Jn^/i 



(0)/,2A0.828 



1 + b, (nf /i^) 



b2 



(All) 



with 



bi = 0.313(n^/i')-°-^^9 1 + 0.607 {Qlh^) 



-0 ,2\ 0.674 



&2 = 0.238 {nl,h^) 



,2\ 0.223 



(A12) 
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The typical value of Zd is about 1021 for f2m^ = 0.276 and h = 0.705 

According to the BAO data from the Two-Degree Field Galaxy Redshift Survey (2dF 
GRS) and the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) 



35|, the distance ratio 



dz at two redshifts z = 0.2 and z = 0.35 is measured to be (i°^g2 = 0.1905 ± 0.0061 and 
d°^Q 35 = 0.1097 ± 0.0036 with the inverse covariance matrix: 



•-^BAO 



30124 -17227 
-17227 86977 



The fo^^ the BAO data is described as 

^2 _ ( th _ obs \ \ ( th 

ABAO ~ i.-'-i.BAO "''i.BAO^ {'^BAO )ij {■^j,BAO 

where Xi,BAo = c^o.ss)- 



X 



obs \ 
j,BAO) ' 



(A13) 



(A14) 



3. Cosmic Microwave Background (CMB) radiation 



The observational data of CMB are sensitive to the distance to the decoupling epoch 
Hence, by using these data we obtain constraints on the model in the high redshift 
regime [z ~ 1000). 



The acoustic scale Ia and the shift parameter TZ 40|] are defined by 

nDA{z^) 



n{z^) 



1 + z^ 



Q'^'Ho (1 + z,) Da{z.) 



where 2;* is the redshift of the decoupling epoch, given by [41 1 



1048 



1 + 0.00124 ffif/i^ 



with 



0.0783 



-0.738 



9l 



-0.238 



1 + 39.5 



0.763 ' 



92 



0.560 



1 + 21.1 (fi[°^/i2 



1.81 



(A15) 
(A16) 

(A17) 

(A18) 



We use the data from Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) obser- 



vations 



|] on CMB. 



The of the CMB data is 



XCMB — (^iJcMB ^i'.CMB) (^CMb)^,- (^jIcMB ^j,CMB) 



(A19) 
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where Xj^cmb = {Ia{z*),TZ{z^), z^) and C^^g is the inverse covariance matrix. The data from 

^ lead to = 302.09, 7^(2*) = 1.725 and = 1091.3 with the 



WMAP7 observations 



inverse covariance matrix: 



'-^CMB 



V 



2.305 29.698 -1.333 
29.698 6825.27 -113.180 
-1.333 -113.180 3.414 



(A20) 



Consequently, the of the observational data is given by 



2 _ ~2 I 2 I 2 
X — XSN + XbAO + XCMB • 



(A21) 



We note that in our fitting procedure, we take the simple method rather than the Markov 
chain Monte Carlo (MCMC) approach such as CosmoMC 42 1. 
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